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of  many  rocket  motors  and  there  has  been  much  concern  In  the  solid  propellant 
rocket  Industry  as  to  the  extent  that  flaws  could  contribute  to  catastroplc  failure. 
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NOMENCLATURE 


A  cross-sectional  area  of  duct 

Br  non-dimensional  blowing  rate 

Cp  constant  pressure  specific  heat 

C  burning  rate  coefficient 

D  propellant  density 

e  roughness  height  for  sand-roughened  pipes 

f  plane  wall  friction  factor 

G  variable  defined  as  2(y-1)T0R/y 

K  variable  to  correct  friction  factor  for  burning 

h  specific  enthalpy  of  gas 

hr  reference  value  of  specific  enthalpy  used  In  a  linear 

approximation  of  the  relationship  between  enthalpy  and 
temperature 
M  Mach  number 

m  mass  flow  rate 

N  exponent  In  burning  rate  law 

P  pressure  at  x 

Pr  arbitrary  reference  pressure 

PQ  stagnation  pressure 

P  non-dimensional  pressure  defined  as  P/Pr 

R  universal  gas  constant 

Re  Reynolds  number  based  on  hydraulic  diameter 

rb  propellant  burning  rate 
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total  surface  area  Sb  +  Sn 
temperature  at  x 
stagnation  temperature 
gas  velocity  at  x 

gas  velocity  at  boundary  layer  of  a  flat  plate 

normal  gas  velocity  at  wall 

distance  along  debond  measured  from  tip 

debond  width  at  x 

ratio  of  specific  heats 

variable  defined  as  (P^  +  G^2)** 

viscosity  of  gas 

variable  directly  proportional  to  mass  flux  and  defined 
as  m/APr 
gas  density 

gas  density  In  boundary  layers  of  flat  plate 
wall  shear  stress 

conditions  at  a  burning  wall 
conditions  at  a  non-burning  wall 
reference  value 
conditions  at  the  wall 
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INTRODUCTION 

Flaws  or  cracks  In  solid  propellant  grains  have  been  responsible 
for  the  failure  of  ftany  rocket  motors.  Thus,  there  has  been  much  con¬ 
cern  as  to  when  and  why  such  flaws  might  propagate  to  the  extent  that 
they  could  contribute  to  catastrophic  failure  (1,2).  Considerable 
work  has  been  conducted  In  an  attempt  to  analyze  the  factors  which 
create  cracks  or  debonds  In  solid  propellants  (3,4).  Among  this  work 
Is  a  method  to  radiographically  detect  flaws  before  the  motor  Is 
fired  (5).  However,  even  the  knowledge  that  the  propellant  grain  Is 
flaw  free  before  firing.  Is  no  assurance  that  flaws  will  not  form 
after  the  motor  Is  fired  due  to  Internal  pressurization. 

Much  of  the  early  work  dealing  with  propellant  cracks  was  carried 
out  to  determine  If  burning  would  propagate  Into  a  crack  existing  In 
a  solid  propellant  grain  (6,7).  The  results  Indicated  definitely 
burning  would  propagate  Into  cracks  or  debonds  and  that  the  rate 
at  which  the  flame  would  propagate  Increased  with  chamber  pressure. 

Before  firing,  defects  In  propellants  are  formed  In  areas  of 
stress  concentration  caused  by  geometry  and  thermal  expansion.  Ther¬ 
mal  effects  arise  when  the  propellant  Is  cast  Into  the  motor  casing 
and  placed  In  an  oven  to  cure  at  temperatures  between  130°F  and  160°F 
depending  on  the  type  of  propellant.  After  the  curing  period  the 
motor  Is  allowed  to  cool  to  room  temperature  and  stresses  are  formed 
due  to  the  difference  In  thermal  expansion  coefficients  of  the  casing 
and  the  propellant.  For  casings  made  of  fiberglass  the  thermal 
stresses  are  somewhat  relieved  due  to  flexure  of  the  casing  whereas 
steel  casings  do  not  yield  as  easily  and  larger  stresses  are  encountered. 
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After  the  motor  Is  fired,  cracks  can  form  due  to  the  developed 
chamber  pressure.  The  increased  chamber  pressure  simply  Increases 
the  stresses  above  those  that  existed  in  the  unfired  state  (1).  Flaws 
formed  before  or  after  Ignition  normally  originate  In  regions  of 
geometric  stress  concentration  such  as  the  tip  of  the  star  pattern. 
Flaws  existing  at  the  tip  of  the  star  pattern  are  extremely  critical 
since  their  propagation  path  to  the  casing  is  minimum.  Tests  have 
been  conducted  at  Thlokol  Chemical  Corporation,  Utah  Division,  to 
determine  the  effects  of  cracks  that  exist  near  the  tip  of  the  star 
pattern  (2).  The  result  Is  usually  propagation  of  the  crack  directly 
outward  to  the  motor  casing  along  the  minimum  length  path. 

On  the  microscopic  level  it  Is  believed  that  cracks  occur  due 
to  an  adhesive  failure  between  the  binder  (or  fuel)  and  the  oxidizing 
agent.  Once  the  adhesion  between  these  particles  and  the  fuel  has 
failed,  small  pockets  are  formed  around  the  oxidizing  agent  particles 
and  stress  concentration  factors  of  approximately  three  are  experienced 
by  the  material.  These  microscopic  flaws  readily  propagate  to  become 
finite  cracks.  Once  a  crack  Is  formed,  It  Is  probable  that  Ignition 
will  occur  within  It  as  soon  as  It  Is  reached  by  the  burning  surface. 
When  this  occurs,  It  Is  of  Importance  to  know  whether  the  pressures 
developed  by  the  combustion  within  the  crack  are  sufficient  to  make 
the  crack  propagate  or  whether  the  flaw  will  merely  add  to  the  ef¬ 
fective  burning  surface  area  of  the  motor. 

In  order  to  study  the  pressures  developed  In  a  burning  flaw. 

It  is  necessary  to  develop  a  qualitative  feeling  for  the  governing 
parameters.  This  includes  a  knowledge  of  what  the  flow  will  look 
like  and  what  parameters  will  affect  it.  If  the  exit  cross-sectional 
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area  of  the  crack  is  small  compared  to  its  burning  area,  the  velocity 
of  the  gaseous  effluent  will  be  much  larger  than  that  external  to 
the  flaw  and  the  flow  will  be  directed  outward  with  no  flow  Into  the 
crack.  If  the  exit  area  Is  larger,  there  may  be  some  partial  Invasion 
of  the  flaw  by  the  external  flow.  Thus,  the  critical  parameters 
which  effect  the  flow  In  the  flaw  and,  hence,  the  pressure  are  the 
propellant  properties,  main  chamber  pressure,  defect  geometry,  main 
chamber  gas  velocity,  propellant  burning  rate,  and  fracture  propa¬ 
gation  velocity.  With  a  knowledge  of  these  governing  parameters, 
one  Is  ready  to  study  the  pressure  distribution  within  the  flaw. 

It  Is  the  main  concern  of  this  study  to  determine  experimentally 
what  the  pressure  distribution  will  be  In  a  propellant  debond  once 
the  burning  has  propagated  Into  It  and  to  compare  the  results  with 
a  one-dimensional  quasi-steady  flow  theory  developed  to  predict 
pressures  In  burning  propellant  flaws.  A  debond  configuration  was 
chosen  to  be  studied  basically  because  of  ease  of  Instrumentation. 

The  results  for  the  debond  configuration  are  not  totally  unrelated 
to  those  for  a  crack,  however.  The  debond  Is  simply  a  crack  with 
propellant  on  one  side  and  the  motor  casing  on  the  other.  Where 
the  crack  has  two  burning  surfaces,  the  debond  has  only  one  which 
makes  pressure  measurements  along  the  non-burning  surface  possible 
without  disturbing  the  flow  pattern.  Results  of  these  measurements 
are  compared  to  the  predicted  values  as  calculated  using  basic  com¬ 
pressible  flow  theory  and  a  computer  model  developed  for  cracks 
or  debonds  In  Reference  8.  Once  the  validity  of  the  predicted  pressure 
values  for  debonds  Is  established,  it  will  be  a  good  assumption  to 
assume  the  theory  holds  true  for  cracks  as  well. 


ONE-DIMENSIONAL  ANALYSIS  OF  A  BURNING  DEBOND 
WITH  FRICTION  AND  MASS  ADDITION 

'Hie  following  analysis  Is  based  on  the  assumption  of  quasi-steady 
state  phenomena.  This  assumption  can  be  justified  on  the  basis  that  the 
propellant  burning  rate  Is  much  slower  than  will  be  the  velocity  of  the 
gaseous  effluent.  Also,  If  mechanical  propagation  of  the  crack  Is  en¬ 
countered,  the  quasi-steady  state  assumption  will  hold  If  the  velocity 
of  mechanical  propagation  Is  much  slower  than  the  gaseous  velocity  or 
If  the  defect  remains  stable  for  a  length  of  time  and  then  propagates 
to  a  new  stable  geometry  which  Is  Instantaneously  Ignited.  Thus,  the 
following  analysis  will  be  to  determine  the  Instantaneous  pressure  dis¬ 
tribution  within  a  burning  propellant  debond. 

The  governing  equations  for  the  state  of  the  effluent  gas  and  the 
basic  geometrical  configuration  are  given  below: 


Continuity 

m  ■  puA 

(1) 

Equation  of  state 

P  -  pRT 

(2) 

Energy 

h  +  u2/2  *  h0 

(3) 

Heat  capacity 

h  -  h  ■  CJT 
r  p 

(4) 
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The  basic  assumptions  Irclude  steady  one-dimensional  adiabatic  flow* 
mass  addition  at  constart  enthalpy  and  negligible  kinetic  energy, 
and  constant  specific  heat  and  viscosity.  The  assumption  of  mass 
addition  at  constant  enthalpy  and  negligible  kinetic  energy  Is  reason 
able  If  the  combustion  zone  Is  assumed  Inflnlt  nimally  thin  and  at 
the  surface.  This  Is  a  good  assumption  for  most  burning  surfaces. 

The  assumption  of  constant  specific  heat  and  viscosity  simply  Im¬ 
plies  that  there  Is  no  chemical  reaction  taking  place  after  the  gas 
Is  added  to  the  stream.  Adiabatic  flow  Is  also  a  good  assumption 
for  solid  propellant  grains  which  release  considerable  energy  and 
are  very  poor  thermal  conductors. 

The  momentum  equation  for  one  dimensional  flow  Is: 


where  A  «  cross-sectional  area;  d$b/dx  •  change  In  surface  area  at 

burning  wall  per  change  In  x;  and  dSp/dx  «  change  In  surface  area  at  a 
non-burning  wall  per  change  In  x.  The  quantity  dy/dx  Is  assumed 

small  to  Insure  the  one  dimensionality  of  the  flow  and  the  x  com- 
ponent  of  momentum  of  the  mass  element,  dm,,  added  over  the  length  dx 
Is  assumed  negligible  due  to  the  assumption  of  normal  burning. 

For  a  debond  the  variation  of  flow  area  and  differential  surface 
area  with  x  are  of  the  form: 

A  ■  y 

dS.  -  dS_  -  dx 
b  n 

The  momentum  equation  contains  two  types  of  wall  shear  stress-shear  at  the 
burning  wall,  Twb,  and  shear  at  the  non-burning  wall,  twr.  Both  types  of  wall 
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shear  stress  can  be  written  In  the  form: 


(7) 


where  the  friction  factor  fnfor  the  non-burning  wall  Is  a  function  of  sur¬ 
face  roughness  and  Reynolds  number.  The  friction  factor  f^  for  a  burning 
wall  Is  dependent  on  Reynolds  number  and  the  mass  addition  rate  which  tends 
to  blow  the  boundary  layer  away  from  the  wall  thus  reducing  the  wall  shear 

stress.  Reference  9  has  reported  a  correlation  between  friction  factor  for 

< 

external  flow  with  pressure  gradient  and  non-dimensional  blowing  rate*  BR, 
where 


B. 


pw  vw  j «r- 

•PC 


(8) 


The  problem  of  flow  In  a  debond  more  closely  resembles  pipe  flow  than  bound¬ 
ary  layer  flow  and  It  would  seem  more  reasonable  to  base  Reynolds  number  on 
hydraulic  diameter  rather  than  distance  from  the  leading  edge.  Doing  this 
the  expression  for  Reynolds  number  becomes: 


Re.ljtt  (9) 

For  one-dimensional  flow  pw  ■  p$  »  p  and  the  normal  velocity  at  the  wall, 
Vw,  can  be  evaluated  by  considering  a  differential  wall  length  dx.  Thus 

dm  -  PVW  dx  (10) 

Combining  Equations  10  and  8  we  obtain 

(11) 

The  earlier  mentioned  correlation  presented  In  Reference  9  may  be  approxi¬ 
mated  to  obtain  the  following  form  for  the  mass  addition  friction  factor  In 
terms  of  the  plane  wall  friction  factor  f: 


where 


VHf 
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H  [s-rra ♦  ^7t]  (,3) 

For  turbulent  flow  a  similar  type  of  correlation  between  fb  and  f  exists 
as  given  In  Reference  10. 

Since  the  case  at  hand  closely  approximates  pipe  flow,  the  fric¬ 
tion  factors  can  be  found  by  the  use  of  a  Moody  diagram  for  rough 
pipes.  The  relationship  for  laminar  flow  Is 

f  ■  |j-  (14) 

For  turbulent  flow  the  friction  factor  for  a  rough  pipe  depends  on  sur¬ 
face  roughness  and  Reynolds  number.  Specific  values  of  f  for  turbulent 
flow  can  vary  between  approximately  0.07  and  0.01  depending  on  the  de¬ 
gree  of  turbulence  and  roughness.  Several  equations  have  been  proposed 
to  characterize  the  behavior  of  f  for  turbulent  flow.  Such  an  equation 
has  been  developed  for  sand  roughened  pipes  In  Reference  11  and  has  the 
following  form: 

1.14  -  0.86  in  §  (15) 

y 

where  e  Is  the  roughness  height  for  sand  roughened- pipes  and  can  be 
compared  to  the  particle  size  of  the  oxidizing  agent  In  a  solid  propellant. 

Substituting  Equations  6,  7,  and  12  Into  the  general  momentum 
equation  one  obtains 

y  H  • fp  1,2  (1r>  *  a?  <■»>  061 

where  f  and  H  are  determined  by  Equations  11,  13,  14,  and  15. 

The  problem  Is  now  defined  In  terms  of  five  basic  equations— energy, 
continuity,  momentum,  heat  capacity,  and  the  equation  of  state.  These 
five  equations  are  expressed  In  terms  of  seven  basic  variables— h,  T,  y,  u. 


s 


P,  p,  and  m.  It  is  now  possible  to  express  any  one  parameter  as  a  func 
tlon  of  any  other  two. 

If  we  use  the  perfect  gas  relation,  Equation  2,  the  property  rela¬ 
tion,  Equation  4,  and  the  continuity  relation,  Equation  1,  to  eliminate 
p,  u,  and  h  from  the  momentum  and  energy  equations,  we  arrive  at 


and 

(Y  -  1)  (Tn2R)  +  (2P2Y)  (Tn2R)  -  2T0Pyt,2  R  •  0  '  (18) 


where 

P  »  P/P  and  n  *  m/P  A 
r  r 

Solution  of  Equation  18  by  the  quadratic  formula  yields 
Tl  9„  _  -P2  ±  Py  (P2  +  62n2)>* 

Tn  r  ^  (19) 

where 


G  «  2(y-1)T0  R/y 

Is  a  constant  parameter  of  the  propellant.  Price  (12)  has  shown1 that  of 

,  \ 

the  two  signs  possible  In  Equation  19  orily  the  plus  sigq  has  physical  sig¬ 
nificance  since  the  minus  sign  corresponds  to  removal  of  mass  In  an  unmix¬ 
ing  process,  which  leads  to  a  decrease  In  entropy.  If  we  define  a  new 
variable;  .  i:  : 

K  •  (P  +  GV)* 

and  combine  Equation >19  with  Equation  17  after  some  rearrangement  we  obtain 


dP  B  ■  Y£  > 

37  y((-V) 


(20) 
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The  burning  rate  of  a  propellant  neglecting  erroslve  burning  may  be 
represented  by  the  empirical  relation 

rb  *  CPN  (21) 

where  rb  Is  the  local  recession  rate  of  the  burning  wall  and  C  and  N 
are  constants  dependent  upon  the  propellant.  The  local  differential 
mass  flux  Is  then  dm  ■  Dr^dx  where  D  Is  the  density  of  the  solid.  The 
local  mass  flux  may  now  be  Integrated  to  obtain 


Drbdx  + 


m. 


(22) 


o 

where  mQ  Indicates  the  presence  of  burning  on  the  head  end  and  has  the 
value 


m„  =  DC  P" 

O  0  0 


(23) 


Combining  Equations  22  and  23  and  the  definition  of  the  variable  n 


m/APr,  one  obtains 

n 

n  ■  r 


ns 


(24) 


where  n0  ■  DCP0N"' 

We  may  now  solve  for  the  pressure  distribution,  or  equivalently  P,  as 
a  function  of  x  by  simultaneously  solving  the  differential  equation  for 
P  (Equation  20)  and  the  Integral  equation  for  n  (Equation  24). 

In  order  to  solve  the  governing  equations,  It  Is  necessary  to  es¬ 
tablish  a  boundary  condition  at  the  flow  exit.  If  the  flow  Is  subsonic 
at  the  exit  plane,  then  the  pressure  must  be  equal  to  the  main  chamber 
pressure,  PCH.  If  the  flow  Is  choked  at  the  exit  plane,  the  boundary 
condition  must  be  Mach  number  equals  one. 

In  terms  of  the  variables  defined  earlier,  the  Mach  number  can  be 
shown  to  be  represented  by  the  relationship 
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H2-ML 

P(Y-l) 


(25) 


A  computer  program  has  been  devised  to  give  the  solution  of  the  govern¬ 
ing  equations  (8).  This  program  will  be  used  to  provide  theoretical  results 
for  comparison  with  the  experimentally  determined  pressure  distributions. 


EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


BASIC  EQUIPMENT 

A  side  view  of  the  model  used  to  experimentally  simulate  a  burning 
debond  with  transducers  In  place  Is  shown  In  Figure  1.  The  dimensions 
of  the  burning  surface  are  three  Inches  wide  by  four  and  one-half  Inches 
long.  The  propellant  slab  Is  0.10  In.  thick  and  can  be  seen  as  a  thin 
layer  on  the  lower  half  of  the  debond  model  In  Figure  1.  These  di¬ 
mensions  were  chosen  basically  to  attempt  to  keep  the  amount  of  pro¬ 
pellant  burned  In  each  test  small  while  keeping  the  dimensions  of  the 
flaw  on  the  same  order  as  those  that  might  exist  In  real  motors.  The 
angle  of  divergence  of  the  debond  model  Is  variable  with  a  maximum 
of  just  over  five  degrees.  This  allows  the  same  model  to  be  used  to 
examine  debonds  with  different  angles  of  divergence.  The  transducers 
shown  In  Figure  1  are  located  at  positions  .19,  1.72,  and  3.72  Inches 
from  the  debond  tip. 

In  order  to  simulate  the  externally  Imposed  chamber  pressure 
tha\  exists  In  a  rocket  motor  experiments  were  conducted  In  a  combus¬ 
tion  chamber  which  Is  shown  In  Figure  2.  The  chamber  Is  approximately 
one  foot  In  diameter  and  three  feet  long  with  a  maximum  working  pressure 
of  1000  psl.  Experiments  were  conducted  with  a  range  of  simulated 
chamber  pressures  using  both  air  and  nitrogen  to  pressurize  the  chamber. 
The  long  narrow  window  on  the  side  of  the  chamber  made  visual  obser¬ 
vation  of  the  tests  possible. 

The  narrow  window  In  the  chamber  and  the  use  of  0.25  Inch  thick 
Plexiglas  side  pieces  In  the  debond  model  allowed  photographic  visu¬ 
alization  of  the  burning  process.  Pyrex  glass  was  Initially  used 
for  the  side  pieces,  however,  It  did  not  work  satisfactorily.  The  high 
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temperature  of  the  burning  propellant  and  the  very  rapid  application 
of  large  amounts  of  heat  caused  the  pyrex  side  pieces  to  crack.  Plex¬ 
iglass  side  pieces  proved  to  work  very  well  In  place  of  the  pyrex.  The 
plexiglass  side  pieces  were  able  to  withstand  the  high  temperature  of 
the  burning  without  melting  away  completely  due  to  the  very  short 
duration  of  the  burn  and  the  ability  of  plexiglass  to  yield  without 
fracture  under  thermal  stress. 

Motion  pictures  were  obtained  at  64  frames  per  second  for  several 
of  the  experiments  using  a  Bolax  16-mm  camera  and  high-speed  Infrared 
film.  Typical  pictures  obtained  from  one  run  are  shown  In  Figure  7. 

The  photos  shown  are  for  the  very  first  portion  of  the  burn  only. 

This  Is  due  to  the  fact  that  smoke  building  up  within  the  chamber 
would  obscure  the  view  of  the  model  after  approximately  0.1  seconds 
In  each  run  and  It  was  Impossible  after  that  to  see  the  burning  through 
the  smoke.  The  photos  that  were  obtained  for  approximately  the  first 
tenth  of  a  second  were  of  acceptable  quality.  With  these  pictures 
one  can  obtain  a  good  qualitative  feeling  for  how  well  the  Ignition 
Is  occurring  and  how  the  propellant  Is  burning  with  time  for  the 
first  part  of  the  burn.  Since  smoke  had  made  visual  recording  of 
how  the  propellant  was  burning  with  time  Impossible*  It  was  necessary 
to  use  those  photos  that  were  obtained  for  the  Initial  portion  of 
the  burn  to  estimate  when  data  could  be  reduced.  It  was  necessary 
to  reduce  the  data  at  a  time  when  the  geometry  had  not  been  changed 
significantly  by  burning. 

As  can  be  seen  In  Figure  2,  tw^  spring-operated  pressure  relief 
valves  are  attached  to  the  combustion  chamber.  These  valves  were 
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used  to  attempt  to  regulate  the  pressure  within  the  chamber.  The 
regulation  achieved  by  the  use  of  the  pressure  relief  valves  was 
relatively  poor  at  low  chamber  pressures  due  to  their  Inability  to 
maintain  large  enough  flow  rates.  Also  the  combustion  gases  that 
were  being  exhausted  through  the  valves  were  highly  corrosive  and 
required  their  frequent  replacement.  Although  the  chamber  pressure 
was  not  held  constant,  the  Increase  In  pressure  that  was  experienced 
by  the  chamber  was  recorded  so  that  the  exact  chamber  pressure  at 
any  Instant  was  known. 

INSTRUMENTATION 

The  pressure  measurements  were  made  at  three  locations  along  the 
debond  and  at  one  location  on  the  pressure  chamber  using  Klstler  Model 
603A  quartz  pressure  transducers.  A  drawing  of  one  of  these  trans¬ 
ducers  Is  shown  In  Figure  3. along  with  a  cross-section  of  the  mounting 
configuration  used.  The  recessed  style  of  mounting  was  chosen  to 
protect  the  transducers  from  the  extreme  heat  of  the  burning  propellant 
(approximately  6000°R).  At  the  high  temperatures  encountered  In  these 
experiments,  the  seal  on  the  transducers  Is  very  Important.  Even  the 
slightest  gas  leaks  can  cause  extremely  high  heat  transfer  rates 
which  may  damage  the  transducer.  A  stainless  steel  sleeve  also  shown 
In  Figure  3  serves  to  provide  the  seal  necessary  to  protect  the  instu- 
ment.  With  this  seal  and  a  snug  fit  the  transducer  Itself  never 
experiences  temperatures  above  Its  working  maximum  of  500°F.  The  1/16- 
Inch  connecting  passage  also  helps  protect  the  transducer  by  removing 
It  from  direct  contact  with  the  hot  flowing  gases  while  at  the  same 
time  allowing  adequate  frequency  response.  With  the  use  of  the  re¬ 
cessed  mounting  configuration  and  given  the  short  duration  of  the 
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(a.)  Transducer  and  Adaptor 
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burn,  these  transducers  will  withstand  gas  temperatures  of  well  above 
6000°F. 

..  After  a  given  number  of  tests  the  transducers  were  each  checked 
against  a  dead  weight  tester  to  determine  if  there  had  been  any  damage 
from  temperature  effects.  The  check  showed  that  the  transducers  had 
not  been  damaged  by  the  tests. 

The  basic  theory  of  operation  of  a  quartz  transducer  begins  with 
pressure  applied  to  the  diaphragm  of  the  transducer  which  Is  converted 
to  a  force  acting  on  the  transducer  crystal.  This  causes  the  trans¬ 
ducer  crystal  to  generate  an  electrical  charge  output  which  Is  pro¬ 
portional  to  the  pressure  Input.  The  basic  sensitivity  of  these 
Instruments  is  "unit  charge  per  unit  pressure"  and  is  expressed  In 
plcocoulombs  per  psl  (pcb/psi).  Because  the  output  of  the  transducer 
Is  a  high- Impedance  charge  signal,  a  charge  amplifier  Is  required. 

The  charge  amplifier  converts  the  high  impedence  charge  signal  from 
the  transducers  to  a  low  Impedence  voltage  or  current  signal  required 
to  drive  the  recording  or  display  equipment. 

The  charge  amplifier  Is  a  dc  voltage  amplifier  with  a  capacitive 
feedback  path  from  the  low-impedence  output  circuit  to  the  hlgh- 
Impedence  Input  circuit.  The  output  voltage  which  results  from  a 
charge  signal  Input  Is  returned  to  the  Input  circuit  through  the 
feedback  capacitor,  In  the  direction  to  maintain  the  Input  circuit 
voltage  at  or  near  zero.  Thus  the  net  charge  from  the  input  circuit 
Is  stored  In  the  feedback  capacitor,  producing  a  potential  difference 
across  it  equal  to  the  value  of  the  charge  divided  by  the  value  of 
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capacitance*  and  this  potential  difference  determines  the  relation¬ 
ship  of  the  output  signal  voltage  to  the  Input  signal  charge.  The 
transfer  characteristic  or  "gain"  of  the  amplifier  depends  only  on 
the  value  of  the  feedback  capacitor  (selected  by  a  range  switch) 
and  the  setting  of  a  precision  voltage  divider  (calibration  factor 
dial),  and  Is  unaffected  by  time*  temperature,  or  line  voltage 
fluctuations  (i3). 

The  charge  signal  from  the  transducer  Is  extremely  sensitive 
to  line  capacitance  or  contaminants  on  the  contact  surface.  For 
this  reason  only  special  low  noise  coaxial  cables  are  suitable  to 
transmit  the  signal  from  the  transducer  to  the  charge  amplifier. 

Any  stray  capacitance  such  as  dirt  or  moisture  In  the  connectors 
will  cause  spurious  signals  due  to  the  breakdown  of  the  transducers' 
high  Insulation  resistance  by  the  contaminant.  The  manufacturer 
suggests  that  freon  solvent,  due  to  Its  excellent  dielectric  pro¬ 
perties,  and  lint-free  paper  wipers  are  the  best  wav  to  keep  the 
connectors  clean.  Other  solvents  such  as  alcohol  tend  to  contami¬ 
nate  Insulators  and  should  not  be  used. 

For  most  applications  the  signal  can  be  fed  directly  from  the 
charge  amplifier  to  the  readout  equipment  without  further  amplifica¬ 
tion.  The  charge  amplifiers  do  have  a  limited  output  capability  and 
In  some  cases  additional  amplifiers  are  required  to  drive  the  readout 
equipment. 

The  equipment  used  for  readout  In  these  experiments  was  a  B  &  F 
model  3006-DL  ultraviolet  oscillograph.  This  Instrument  utilizes 
light  sensitive  paper  and  a  zenon  lamp  to  produce  a  trace.  The 
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oscillograph  has  a  capacity  of  twelve  channels*  each  channel  rqoulr- 
Ing  a  galvanometer  for  operation.  Only  four  channels  were  required 
for  this  work,  however. 

The  galvanometers  that  were  utilized  here  consist  of  simply 
a  small  mirror  suspended  by  a  coll  of  very  fine  wire.  The  deflection 
In  the  trace  Is  achieved  by  the  magnetic  field  applied  to  the  coll 
when  current  Is  passed  through  It  which  subsequently  causes  rotation 
of  the  mirror  In  proportion  to  the  magnitude  of  the  applied  magnetic 
field.  It  was  found  that  the  current  output  of  the  charge  amplifiers 
Is  somewhat  limited  and  was  Insufficient  to  drive  the  galvanometers 
full  scale  at  the  voltage  levels  encountered.  Thus  it  was  necessary 
to  place  an  operational  D.  C.  amplifier  In  line  between  each  of  the 
charge  amplifiers  and  the  oscillograph  to  provide  the  necessary 
current.  For  convenience,  the  amplifiers  were  chosen  with  a  gain 
setting  of  unity  In  order  that  their  current  output  capabilities  could 
be  taken  advantage  of  without  amplifying  any  noise  that  might  exist 
In  the  system.  A  schematic  of  the  Instrumentation  circuit  for  one 
channel  Is  shown  In  Figure  4. 

An  ultra-violet  recorder  was  chosen  because  the  extreme  transient 
nature  of  the  tests  required  that  the  output  device  have  a  very  fast 
writing  speed  and  also  a  very  high  frequency  response.  The  writing 
speed  of  the  light  sensitive  recorder  used  Is  on  the  order  of  a  max¬ 
imum  of  60,000  Inches  per  second  and  the  frequency  response  of  the 
galvanometers  that  weVe  chosen  was  at  least  1000  Hertz.  This  com¬ 
bination  of  fast  writing  ability  and  quick  response  was  more  than 
adequate  to  record  the  events  accurately  without  any  appreciable  lag. 
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In  order  that  the  Ignition  process  be  recorded  correctly,  It 
was  required  that  the  paper  drive  on  the  recorder  be  started  just 
slightly  before  the  Ignition  took  place.  This  would  be  no  problem 
to  achieve  manually  for  slow  paper  speeds,  but  at  the  faster  speeds 
used  It  would  mean  wasting  large  amounts  of  light  sensitive  paper 
which  was  quite  expensive.  To  save  on  paper  and  make  the  Ignition 
simple  a  time-delay  relay  was  Incorporated  In  the  Ignition  circuit. 

A  single  throw  double-pole  switch  was  used  to  start  the  paper  drive 
and  simultaneously  energize  a  timed  relay  which  was  set  to  provide 
current  to  the  Ignition  wire  and  an  event  marker  on  the  oscillograph 
after  approximately  a  two  second  delay.  This  circuit  made  It  possible 
to  begin  the  entire  event  by  simply  throwing  a  single  switch.  The 
Ignition  circuit  Is  shown  In  block  diagram  form  In  Figure  5a  and 
the  wiring  diagram  Is  shown  In  Figure  5  b. 

The  paper  drive  was  set  at  4  Inches  per  second  and  the  paper 
was  marked  with  vertical  timing  lines  at  0.1  second  Intervals.  The 
event  marker  made  It  possible  to  determine  the  time  elapsed  between 
Ignition  and  the  Initial  displacement  of  the  recorder  trace. 

IGNITION  OF  A  SOLID  PROPELLANT  SURFACE 

As  has  been  mentioned  the  surface  of  the  propellant  In  the  debond 
model  was  three  Inches  wide  and  four  and  one-half  Inches  long.  In 
order  that  the  tests  be  meaningful  and  comparable  with  the  computer 
analysis.  It  was  required  that  the  entire  surface  be  burning  at  the 
time  of  the  pressure  measurement.  This  made  It  necessary  to  Ignite 
the  propellant  over  the  entire  surface  "Instantaneously."  The  word 
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Instantaneously  being  used  to  mean  that  the  entire  surface  should  be 
burning  before  any  appreciable  burning  normal  to  the  surface  has  oc¬ 
curred.  The  propagation  rate  of  burning  over  a  free  surface  Is  much 
faster  than  the  burning  rate  normal  to  the  surface  and  this  effect 
Is  Increased  noticeably  at  higher  pressures.  Several  methods  were 
Investigated  to  attempt  to  achieve  suitable  Ignition  over  the  entire 
surface. 

The  goal  was  to  achieve  very  rapid  Ignition  of  the  surface  at 
atmospheric  pressure  so  that  the  Ignition  under  Increased  pressures 
could  be  safely  assumed  to  approach  that  of  Instantaneous  Ignition 
over  the  entire  surface. 

Among  the  Methods  of  achieving  Instantaneous  Ignition  of  a  sur¬ 
face  Is  that  of  painting  the  surface  of  the  propellant  with  a  mixture 
of  red  phosphorus  powder  and  an  organic  binder  such  as  common  Duco 
cement.  This  method  was  among  the  first  tried  and  did  not  give  satis¬ 
factory  results.  The  reason  being  perhaps  that  the  red  phosphorus 
which  had  been  used  had  been  exposed  to  moisture  and  had  lost  Its 
flammability  characteristics.  The  next  method  attempted  was  simply 
roughening  the  surface.  While  this  did  tend  to  make  the  flame  spread 
more  rapidly  across  the  surface,  the  propagation  rate  of  the  burning 
was  still  not  rapid  enough  to  fulfill  the  requirement  of  Instantaneous 
Ignition. 

Another  common  method  of  achieving  rapid  surface  Ignition  Is 
spreading  a  highly  flammable  powder  over  the  surface.  The  first  sub¬ 
stance  that  was  Investigated  was  simply  powdered  propellant.  The  pro¬ 
pellant  powder  Ignited  extremely  fast  but  produced  only  enough  heat 
to  Ignite  the  solid  body  of  propellant  In  a  few  locations.  The  flash 
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of  the  propellant  powder  was  on  the  order  of  the  Ignition  speed  that 
was  being  sought  although  It  did  not  Ignite  the  propellant  properly. 

The  problem  was  now  to  find  a  powder  that  would  flash  when  Ignited  and 
produce  enough  heat  to  start  the  propellant  burning  everywhere. 

A  method  to  retain  the  heat  produced  by  a  powder  that  burns 
extremely  fast  Is  mixing  Iron  filing  or  aluminum  paint  flake  with  the 
powder.  Iron  filings  were  mixed  with  both  propellant  powder  and  black 
gun  powder  and  tested  for  Ignition.  Neither  method  seemed  to  work 
although  the  black  powder  gave  better  results  than  any  of  the  method 
previously  attempted.  The  Iron  filings  used  were  quite  coarse  and 
perhaps  If  they  were  smaller,  the  test  would  have  given  better  results. 

It  was  finally  discovered  that  metals  such  as  titanium,  zerconlum, 
and  magnesium  when  ground  to  very  fine  powders  and  mixed  with  an  ox¬ 
idizer  will  burn  extremely  fast  and  give  off  large  amounts  of  heat  at 
the  same  time.  The  powdered  metals  must  have  particle  sizes  on  the 
order  of  thirty  to  fifty  microns  or  smaller.  The  smaller  the  particles, 
the  faster  will  be  the  burning  rate  and  the  smaller  the  time  for  heat 
transfer  to  the  propellant.  Of  the  many  mixtures  possible,  It  was 
found  that  a  mixture  of  seventy  per  cent  titanium  and  thirty  per  cent 
ammonium  perchlorate  powder  worked  very  well.  The  titanium  used  had 
an  average  particle  size  of  five  microns  and  the  ammonium  perchlorate 
particles  were  fifteen  microns  on  the  average.  Larger  particles  up 
to  twenty  or  thirty  microns  average  could  be  used  with  the  same  results. 

The  Ignition  procedure  was  then  defined  as  follows.  A  very  fine 
nlchrome  wire  was  attached  to  the  surface  of  the  propellant  In  a 
pattern  as  shown  In  Figure  6.  On  top  of  the  wire  and  the  propellant 
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there  was  then  spread  a  very  fine  layer  of  the  aforementioned  mixture 
of  titanium  Ignition  powder.  It  was  very  Important  that  the  powder 
be  spread  evenly  over  the  entire  surface  so  that  the  heat  from  the 
flash  would  reach  every  section  of  the  surface  evenly.  With  the  wire 
and  powder  In  place  the  debond  model  was  ready  for  Ignition. 

Several  trial  tests  were  run  and  motion  pictures  taken  to  try  to 
determine  If  the  Ignition  was  Indeed  as  good  as  It  appeared  to  be. 
Figure  7  shows  three  frames  blown  up  from  these  pictures.  The  top 
picture  occurred  0.0625  seconds  afte**  Ignition  contact  and  the  other 
two  followed  In  sequences  of  1/64  of  a  second.  The  flow  pattern 
developed  In  the  top  picture  was  the  first  evidence  of  burning  to 
occur  In  the  sequence.  In  the  bottom  picture  the  propellant  slab 
can  be  seen  with  little  sign  of  geometry  change  due  to  burning.  The 
tip  area  Is  obscured  by  smoke  and  the  ring  section  of  the  model  In 
some  of  the  photographs.  Although  the  pictures  taken  were  not  con- 
elusive,  It  seems  that  the  propellant  was  being  Ignited  quite  evenly 
and  very  rapidly.  From  visual  observation  of  the  Ignition  phenomena 
of  several  tests.  It  was  concluded  that  the  method  worked  very  well 
In  getting  the  propellant  Ignited  over  the  entire  surface  evenly  and 
rapidly  and  also  the  Ignition  was  very  repeatable  from  test  to  test. 


DISCUSSION  OF  RESULTS 


The  experimental  results  are  shorn  In  Figures  8  through  14.  Com¬ 
parison  with  numerical  solutions  are  Indicated.  Experiments  were  con¬ 
ducted  for  debond  angles  of  4.75  degrees  and  3.17  degrees  and  1.90  degrees 
and  the  data  reduced  at  chamber  pressure  of  170,  245,  300,  332,  400,  550, 
and  600  p$1.  The  data  Is  plotted  with  an  error  band  of  +  ten  per  cent 
to  give  the  reader  an  Idea  of  the  magnitude  of  the  experimental  error 
Incurred  In  each  case.  The  details  of  the  geometries  are  given  for  each 
case  along  with  the  plot  of  the  results. 

The  tip  width  specified  in  each  case  Is  a  necessary  geometrical 
consideration  In  order  to  eliminate  the  singularity  that  exists  In  the 
mathematical  analysis  for  a  tip  width  of  zero.  As  can  be  seen 
trend  of  the  curves  predicted  by  the  theory*  the  pressure  approaches 
Infinity  asymptotically  for  small  values  of  tip  width  and  x  close  to 
zero.  Since  the  pressure  cannot  logically  be  Infinite  at  the  tip  and 
also  to  make  a  solution  possible  mathematically,  It  becomes  necessary  to 
assume  a  tip  width  other  than  zero.  This  assumption  Is  valid  on  the 
basis  that  a  truly  sharp  crack  does  not  exist  In  nature.  Determination 
of  the  tip  width  to  be  used  Is  based  somewhat  on  a  knowledge  of  the  con¬ 
figuration  to  be  studied  and  an  analysis  of  how  the  tip  width  will  effect 
the  theoretical  pressure  distribution  within  the  debond. 

The  effect  of  a  change  In  tip  width  on  the  predicted  values  of  tip 
pressure  and  also  pressure  distribution  Is  shown  In  Figure  15  for  a  typical 
debond  configuration.  As  can  be  seen  from  the  figure,  the  distribution 
of  pressure  within  the  flaw  Is  little  affected  for  varying  tip  widths 
below  0.001  of  an  Inch  and  beyond  a  point  approximately  0.6  Inches  away 
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from  the  tip.  For  points  closer  to  the  tip*  the  pressure  will  be  some¬ 
what  affected  by  the  tip  width  assumption.  It  can  be  seen  from  the 
figure  that  the  changes  Induced  In  the  debond  pressure  distribution 
by  making  the  tip  width  smaller  are  limited  to  the  Immediate  vicinity 
of  the  tip  and  that  they  approach  the  case  where  only  the  tip  pressure 
Itself  Is  affected  by  a  change  In  tip  width.  With  this  In  mind  and  with 
some  consideration  of  the  Increased  computer  time  required  for  smaller 
tip  widths,  a  tip  width  of  0.001  Inches  was  chosen  to  be  used  In  the 
theoretical  predictions.  Once  the  tip  width  was  chosen.  It  was 
necessary  to  define  the  propellant  properties. 

The  propellant  utilized  In  these  experiments  Is  a  double-based  pro¬ 
pellant  loaded  with  aluminum.  The  parameters  that  are  required  for  the 
analysis  Include  burning  temperature,  burning  rate  coefficient,  burning 
rate  exponent,  viscosity  of  gaseous  effluent,  propellant  density,  and 
the  ratio  of  specific  heats  for  the  gaseous  effluent  (1,  2). 

The  data  was  reduced  for  each  experiment  by  simply  picking  the  chamber 
pressure  for  the  case  that  was  desired  to  be  studied  and  then  finding 
at  which  point  the  recorded  chamber  pressure  trace  crossed  this  value  as 
long  as  the  desired  chamber  pressure  was  reached  In  the  first  half 
second  of  the  test.  This  latter  restriction  was  placed  on  the  experi¬ 
ment  due  to  the  fact  that  thermally  Induced  drift  exists  In  the  pressure 
transducers.  This  technique  of  obtaining  data  at  a  specified  pressure 
was  also  necessary  due  to  the  continuously  changing  chamber  pressure.  Values 
of  pressure  were  read  from  the  other  three  recorded  traces  at  points 
corresponding  to  that  of  the  desired  chamber  pressure.  The  one-half 
second  limitation  was  also  necessary  to  pick  a  chamber  pressure  In  order 
to  be  sure  that  the  burning  had  not  changed  the  geometry  appreciably. 
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Also,  readings  could  not  be  taken  too  soon  after  the  Ignition  switch 
had  been  thrown  In  order  that  stable  burning  might  have  time  to 
establish  Itself.  This  limited  the  data  that  could  be  meaningfully 
reduced  to  that  which  occurred  approximately  between  0.20  and  0.5 
seconds  after  Ignition.  The  ayerage  time  between  electrical  Ignition 
contact  and  Initial  trace  displacement  was  approximately  0.05  to  0.12 
seconds.  This  delay  was  most  likely  due  to  heating  of  the  Ignition 
wire  and  also  the  time  required  for  flow  to  begin  moving  outward  from 
the  debond  after  the  surface  had  been  Ignited  (see  Figure  7).  The 
delay  encountered  was  very  repeatable  from  one  experiment  to  the  next 
and  lent  Itself  to  the  accuracy  of  the  Instantaneous  Ignition  assumption. 

The  data  represented  In  Figures  10  and  11  was  taken  from  the  same 
four  experiments.  The  data  In  Figure  10  for  a  chamber  pressure  of  550  psl 
was  In  the  Ideal  range  for  data  reduction  as  explained  In  the  previous 
paragraph  whereas  the  data  In  Figure  11  for  a  chamber  pressure  of  600  psl 
was  taken  at  a  later  point  after  the  chamber  pressure  had  risen  an  add¬ 
itional  50  psl.  A  comparison  of  the  error  encountered  In  each  case  will 
yield  larger  errors  for  the  case  taken  at  the  chamber  pressure  of  600  psl 
as  compared  to  the  one  at  550  psl.  This  tells  one  that  the  error  Is 
Increasing  as  the  burning  continues  which  Is  what  would  be  expected  due 
to  the  change  In  geometry.  On  the  other  hand,  taking  data  too  soon  after 
Ignition  will  yield  even  larger  errors  due  to  the  uneven  burning  that  occurs 
In  the  first  0.1  to  0.2  seconds  after  electrical  Ignition  contact. 

The  data  In  Figures  8  and  9  were  reduced  from  separate  sets  of  ex¬ 
periments  with  each  set  of  data  being  taken  approximately  between  0.20 
and  0.3  seconds  after  Ignition.  The  time  required  for  the  desired  chamber 
pressure  to  be  reached  In  each  case  varied  due  to  Inconsistencies  In 


39 


the  operation  of  the  mechanical  pressure  relief  valves. 

As  can  be  seen  on  any  one  of  the  graphs  of  experimental  results,  the 
number  of  data  points  Is  greatest  near  the  exit,  smaller  near  the  center, 
and  even  smaller  near  the  tip.  The  reason  for  this  being  that  In  some 
particular  cases  one  of  the  channels  would  give  either  no  output  at  all 
or  would  overscale  completely.  A  galvanometer  failure  and  this 
problem  were  responsible  for  the  small  number  of  points  presented  In  each 
of  Figures  12,  13,  and  14.  The  channels  that  did  function  properly 
however,  show  very  good  agreement  with  the  theory  In  every  case. 

An  explanation  of  this  behavior  would  be  contaminants  on  the  connectors 
or  transducers  themselves.  With  the  extreme  sensitivity  to  contaminants 
of  these  Instruments  as  explained  earlier  the  slightest  bit  of  dirt  or 
moisture  would  render  the  output  meaningless.  As  for  the  Increased 
frequency  of  bad  data  for  the  tip  transducer,  the  explanation  would  be 
that  the  effect  of  contaminants  must  be  a  function  of  the  charge  output 
of  the  transducer  or  equivalently  the  magnitude  of  the  pressure  being 
measured.  This  seems  reasonable  on  the  basis  that  the  contaminants  act 
as  stray  capacitance  In  the  line  with  the  effect  of  this  stray  capacitance 
being  dependent  on  the  output  of  the  transducer.  Hence,  the  larger  the 
output  of  the  transducer,  the  larger  the  chance  of  the  output  being 
affected  by  contaminants.  Also  It  was  found  later  on  In  the  experimental 
work  that  there  were  a  few  bad  connections  existing  In  the  circuitry. 

These  bad  connections  were  of  the  type  that  would  either  work  properly 
or  not  at  all.  Hence,  there  existed  a  possibility  of  a  particular 
channel  not  working  at  all  for  any  given  test  due  to  a  loss  of  electrical 
contact.  In  some  cases  the  contacts  were  good  and  a  particular  channel 
would  fall  to  give  a  meaningful  output  due  to  other  effects. 


ERROR  ANALYSIS 


The  average  errors  Incurred  In  the  experimental  work  varied  from  a 
low  of  approximately  0.1  per  cent  to  a  high  of  19.8  per  cent  where  aver¬ 
age  per  cent  error  Is  taken  to  mean  the  difference  between  the  theore¬ 
tical  prediction  and  the  ayerage  of  the  experimental  values  divided  by 
the  magnitude  of  the  theoretical  prediction.  A  summary  of  the  average 
errors  encountered  for  some  of  the  cases  studied  Is  given  below. 


Per  Cent  Average  Error 


Chamber 

Pressure, 

Psl 

m 

Center 

Exit 

170 

5.5 % 

4.2% 

0.1% 

300 

14.4% 

10.8% 

0.1% 

550 

11.4% 

2.7% 

2.1% 

600 

19.8% 

2.2% 

3.5% 

Where  the  averages  as  shown  above  reflect  reasonable  errors  In  most 
cases,  a  considerable  amount  of  scatter  can  be  seen  In  the  data  especially 
In  Figures  8  and  9  and  the  scatter  Is  seen  In  Increase  for  points  close  to 
the  tip.  The  explanation  for  this  behavior  would  most  likely  lie  In  the 
transducer  and  cables.  As  explained  before  contaminants  could  be  res¬ 
ponsible  for  this  error. 

The  experimental  procedure  called  for  extreme  care  In  attempting  to 
keap  the  transducers  and  connectors  clean.  Freon  T-F  spray  was  used  to 
clean  the  contact  surfaces  between  each  experiment  and  the  Inside  of  the 
pressure  chamber  was  washed  frequently.  It  Is  believed,  however,  that 
due  to  the  large  amount  of  contaminants  produced  by  the  burning,  there  Is 
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little  possibility  of  excluding  these  particles  of  dust  or  moisture  from 
finding  their  way  Into  the  electrical  contacts.  The  real  problem  arises 
after  the  burn  Is  made  and  the  connections  between  the  low  noise  cable 
and  the  transducers  are  broken  In  order  to  remove  the  debond  model  from 
the  bomb  to  be  cleaned  and  recharged  with  new  propellant.  At  this  time 
the  Inside  of  the  chamber  Iscompletely  covered  with  combustion  products 
from  the  burn  Including  some  moisture.  With  cautious  preparation  the 
effects  of  contaminants  can  be  reduced  somewhat  but  It  Is  not  reasonable, 
considering  the  large  amount  of  material  produced  by  burning,  that  their 
effect  could  be  removed  100  per  cent. 

Contaminants  could  also  have  been  responsible  for  errors  In  one 
other  way.  By  Inspection  of  the  debond  model  after  a  burn  It  was  noticed 
that  deposits  of  a  slag-like  material  had  built  up  In  small  mounds  on  both 
sides  of  the  debond,,  see  Figure  16  and  Figure  17.  While  these  deposits 
present  no  problem  to  the  flow  nor  will  they  affect  the  pressure  distribution 
developed  within  the  model,  they  do  present  a  possible  source  of  error  In 
the  readout.  Although  It  cannot  be  substantiated.  It  seems  possible  that 
on  certain  runs  the  1/16  Inch  connecting  passages  between  the  burning 
propellant  and  the  transducers  could  be  clogged  or  partially  blocked  by 
the  buildup  of  this  slag.  This  blocked  passage  would  then  be  unable  to 
transmit  the  pressure  correctly  and  the  output  of  the  transducer  would  be 
In  error.  This  then  becomes  another  possible  explanation  for  erroneous  or 
Incorrect  data. 

Experimental  errors  could  have  also  developed  from  any  one  of  several 
other  factors  along  with  those  just  discussed.  For  example,  the  burning 
rate  law  which  Is  used  In  the  mathematical  analysis  could  be  questioned 
somewhat.  The  burning  characteristics  of  a  solid  propellant,  particularly 
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TP-H1011,  vary  considerably  with  pressure.  The  burning  rate  law  used 
here  was  arrived  at  experimentally  by  Thlokol  engineers  (1,  2)  at  a 
chamber  pressure  of  about  724  psl.  There  Is  some  question  as  to  the 
exactness  of  this  relation  used  for  lower  pressure  ranges.  Since  It  Is 
believed  that  propellant  of  this  type  will  burn  with  more  stability  at 
higher  pressures,  the  burning  rate  law  then  becomes  a  possible  source 
of  discrepancy  when  used  at  lower  pressures. 

Along  this  same  line, error  could  be  Induced  Into  the  analysis  by 
the  specification  of  propellant  properties  such  as  burning  temperature, 
density,  viscosity  of  effluent  gases,  and  the  ratio  of  specific  heats. 
These  properties  vary  generally  from  mix  to  mix  and  only  average  values 
can  be  specified.  Where  the  burning  temperature  Is  more  a  function  of  the 
burning  conditions,  particularly  pressure,  the  other  properties  can  vary 
with  age  and  condition  of  the  propellant  specimen.  The  analysis  assumes 
that  the  temperature  of  the  burning  1$  constant  where  In  reality  It  Is  a 
weak  function  of  pressure.  Knowing  the  correct  burning  temperature  Is  a 
very  Important  factor  In  the  analysis. 

That  the  burning  Is  more  stable  at  higher  pressures  can  be  seen  In 
the  reduction  of  scatter  In  the  experimental  data  at  the  higher  pressures 
compared  to  those  at  the  lower  chamber pressures.  Also,  since  the  burning 
rate  Increases  with  pressure,  It  Is  more  likely  that  the  Ignition  will 
occur  In  a  more  stable  fashion  at  higher  pressures.  Countering  the  favor* 
able  effects  of  higher  pressures  Is  the  fact  that  the  Increased  burning 
rate  makes  the  problem  more  transient  In  nature  and  the  geometry  more 
difficult  to  define  after  the  burning  has  begun. 

An  additional  difficulty  In  defining  the  Initial  geometry  of  the 
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debond  Is  the  variation  encountered  In  the  physical  thickness  of  the 
propellant  slab.  The  machining  process  used  to  cut  the  slabs  could 
only  guarantee  a  tolerance  In  thickness  of  +.0.03  Inches  due  to  the" 
viscoelastic  behavior  of  the  propellant  while  being  cut.  The  effect 
of  this  variation  In  thickness  along  with  the  thickness  of  the  layer 
of  adhesive  used  to  bond  the  propellant  to  the  model  beomces  signifi¬ 
cant.  The  error  Induced  by  this  variance  In  thickness  Increases  as 
the  angle  of  divergence  of  the  debond  decreases  and  errors  of  as  high 
as  ten  per  cent  In  the  specification  of  geometry  can  be  Incurred  due  to 
this  effect  alone.  The  effect  also  will  have  larger  magnitudes  near  the 
tip  where  a  variance  of  geometry  on  the  order  of  +  0.03  Inches  can  make 
as  high  as  100  per  cent  change  In  the  debond  width. 

In  Figure  16  a  dark  area  can  be  seen  which  Is  approximately  one 
Inch  wide  and  at  the  end  of  the  model  where  the  tip  exists.  This  burn 
was  made  at  an  angle  of  3.17  degrees  and  the  dark  area  In  question  was 

a  result  of  the  propellant  slab  being  thicker  than  0.1  Inches  which 
caused  the  propellant  co  be  In  contact  with  the  top  surface  of  the  debond 

completely  covering  the  passage  to  the  tip  transducer.  This  caused  a  de¬ 
lay  In  the  tip  transducer's  response  for  a  time  sufficient  to  allow  burning  to 
open  up  the  crack  and  uncover  the  passage  to  the  transducer.  This  drastic 
change  In  the  geometry  made  this  test  less  meaningful  and  the  tip  pressure 
totally  meaningless.  This  effect  was  present  to  a  lesser  degree  In  other  runs. 

Two  other  possible  sources  of  error  were  Investigated  to  determine 
their  contribution.  These  were  possible  temperature  effects  on  the 
transducers  and  also  the  possibility  of  the  effect  of  Ionized  gasses 
near  the  transducer.  In  order  to  determine  If  changes  In  the  temperature 
of  the  transducer  Itself  would  effect  the  output*  an  acetlyene  torch 


was  used  to  simulate  the  hot  flowing  gasses  that  exist  In  the  debond 
while  at  the  same  time  creating  no  appreciable  pressures  on  the  trans¬ 
ducers.  This  test  was  made  with  the  lower  half  of  the  debond  model 
taken  off  so  that  the  hot  gasses  could  flow  freely  over  the  surface  where 
the  transducers  were  mounted.  The  test  showed  that  there  was  little  or 
no  effect  due  to  the  existence  of  the  hot  gas  but  that  there  was  some 
appreciable  drift  due  to  conduction  through  the  plate  to  the  transducer. 
This  drift  was  on  the  order  of  an  equivalence  of  SO  psl  but  did  not  occur 
until  approximately  two  to  three  seconds  after  application  of  the  hot 
torch.  However,  the  6000°F  temperature  of  combustion  of  the  p^pellant 
within  the  debond  could  cause  a  larger  drift.  Thus,  it  was  concluded 
that  temperature  effects  due  to  conduction  thru  the  plate  although  present 
In  significant  magnitude  were  not  present  during  the  critical  test  period 
of  0.5  seconds  and  therefore  had  no  bearing  on  the  test  results. 

To  test  for  the  effects  of  Ionized  gases  the  same  test  with  the 
torch  was  conducted  with  the  recessed  cavity  and  connecting  passage 
filled  with  bearing  grease.  This  way  no  Ionized  gasses  could  come  In 
contact  with  or  get  near  the  transducers.  The  results  of  the  test  were 
the  same  as  those  obtained  with  no  grease  In  the  recessed  cavity  and 
connection  passage.  Thus  It  could  be  concluded  that  Ionized  gasses  had 
no  effect  on  the  output  of  the  transducers. 

If  one  studies  the  trend  of  the  experimental  results  compared  to 
the  theory  for  the  tip  transducer  It  will  become  evident  that  experimental 
values  lie  above  predicted  values  In  most  of  the  cases.  If  the  assumed 
tip  width  were  changed  from  .001  Inches  to  .0001  Inches  a  change  In  the 
predicted  pressures  near  the  tip  would  occur  as  seen  In  Figure  1  .  This 
Increase  In  the  magnitude  of  the  predicted  values  near  the  tip  would 
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significantly  reduce  the  variance  between  experimental  and  theoretical 
values.  Thus  the  values  of  average  error  for  the  tip  transducer  as 
given  earlier  In  this  section  are  high  and  would  be  reduced  significantly 
by  use  of  a  smaller  assumed  tip  width.  The  data  Is  left  In  terms  of  a 
tip  width  of  0.001  of  an  Inch  due  to  difficulty  In  the  computer  Iteration 
scheme  encountered  when  the  tip  Is  assumed  smaller. 

With  some  consideration  of  the  large  number  of  factors  which  con¬ 
tribute  to  the  possibility  of  errors  In  the  experimental  data,  the  de¬ 
viations  that  were  encountered  become  of  reasonable  magnitude.  Even  the 
large  amount  of  scatter  that  was  seen  In  some  cases  seems  reasonable 
given  the  factors  which  could  cause  behavior  of  this  typo.  With  all  of 
the  errors  mentioned,  there  has  as  of  yet  been  no  mention  of  assumptions 
given  In  the  theoretical  analysis  such  as  one-dimensional  flow,  etc. 

These  assumptions  could  also  make  induced  error.  This  consideration  along 
wltt  the  experimental  error  analysis  makes  the  results  feasible  on  the 
basis  of  an  engineering  error  analysis. 

i 

i 
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I 
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CONCLUSIONS 


The  experimental  work  conducted  gives  good  substantiation  of  the 
pressure  distribution  In  a  debond  as  calculated  using  one-dimensional 
compressible  flow  theory.  The  cases  studied  cover  a  range  of  chamber 
pressures  from  170  psl  to  600  psl  and  a  range  of  debond  angles  from  1.90 
degrees  of  4.75  degrees.  The  range  of  angles  studied  cover  completely 
those  angles  allowed  under  the  one-dimensional  flow  assumption  and  the 
chamber  pressures  cover  the  range  of  pressures  allowed  for  safety  by  the 
experimental  equipment. 

It  would  be  reasonable  to  conclude  at  this  time  that  the  analytical 
and  numerical  methods  used  here  to  predict  pressure  distributions  In 
burning  propellant  debonds  and  cracks  are  valid  for  angles  below 
approximately  five  degrees  and  chamber  pressures  In  the  range  studied. 
Also  there  1$  good  reason  to  believe  that  the  results  would  compare  more 
favorably  for  tests  run  at  higher  chamber  pressures. 

With  a  working  method  to  predict  pressures  In  flaws  It  Is  now 
possible  to  couple  one-dimensional  flow  analyses  with  a  study  of  the 
mechanical  properties  of  propellants  and  predict  when  failure  will 
occur  due  to  an  existing  flaw  In  a  solid  propellant  grain. 
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APPENDIX  A 


A  listing  of  the  computer  program  used  to  solve  equations  20 
and  24  numerically  Is  given  In  this  appendix.  This  program  Is  a 
slightly  revised  version  of  the  program  developed  In  Reference  8. 

The  baste  solution  technique  Is  the  same  with  changes  occurring  In 
some  areas  to  make  the  program  better  fit  the  debond  configuration 
being  studied. 

The  program  basically  consists  of  two  loops  with  the  solution 
technique  being  an  Iterative  procedure.  The  basic  method  which 
the  program  utilizes  Is  that  of  Initially  quesslng  a  tip  pressure 
and  then  based  on  that  guess  the  distribution  of  pressure  within  the 
entire  debond  Is  calculated.  Once  this  Is  done  the  exit  boundary 
conditions  can  be  checked  and  a  new  value  for  tip  pressure  can  be 
calculated  based  on  the  old  value  and  the  error  encountered  In  the 
boundary  condition  check.  The  procedure  Is  repeated  until  the  boun¬ 
dary  condition  check  and  the  proper  distribution  of  pressure  Is 
found. 

The  basic  object  of  having  two  such  similar  loops  Is  that  loop  1* 

In  general,  will  differentiate  the  choked  from  the  non-choked  flows 
or  In  other  words  loop  1  attempts  to  match  the  boundary  condition 
that  P(N)  *  PC  and  V(N)  «  1.  The  second  loop  basically  takes  care 
of  choked  flows  l.e.  V(N)  «  1. 

Table  A-l  gives  a  comparison  of  the  variables  used  In  the  computer 
program  and  the  analyses.  For  a  more  detailed  explanation  of  the  program, 
a  flow  chart,  and  a  listing  that  will  handle  cracks  as  well  as  debonds 
the  reader  should  see  Appendix  B  In  Reference  8. 
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TABLE  A1.  Nomenclature  and  Units  Used  In  Computer  Program 


Analysis 

Comouter 

Units 

Br 

BR(I) 

none 

C 

C 

ft2n+1 /sec-1 bfn 

0 

0 

1bm/ft3 

f 

FC 

none 

G2 

F 

(sec-lbf/lbm)2 

H 

HBR 

none 

N 

V(I) 

none 

n 

XP 

none 

P 

P(I) 

psl 

Po 

PO 

psl 

Pch 

PC 

psl 

V 

Z(I) 

none 

F' 

ZP 

1/1n. 

R 

R 

ft-lbf/1bm°R 

Re 

RE(I) 

none 

To 

TO 

°R 

X 

S(I) 

In. 

AX 

H 

In. 

y 

Y(I) 

In. 

y* 

YP(I) 

none 

y 

GM 

none 

c 

X{I) 

none 

U 

1bf/sec/ft2 

n 

T(I) 

sec/ft 

% 

THO 

none 
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DIMENSION  Y (2000) r Yp(2000) rS(2000)  »2<2000> »X<2000> , 

♦  V (2000) 

DIMENSION  P(2000) »T(2000) »RE<200Q) »BR(2000) » JJ<3) 
real  ll 

b  READ  3#  YN»Lu»PC» YO 

DATA  DlF»ERROR#H#L»MAX/.Ol»i.».Ql»10»lOO/ 

N  -  451 

na=n 

c**********************************************************c 

c  LlsDEdOND  length 
C  YNZDEBOND  exit  width 
C  PCsCHAMtiER  pressure 
C  YOSTIP  WIDTH 

c  computer  parameters 

c  unallowable  error  between  actual  mach  no  and  mach  1,0 
C  tHRoRsERROR  IN  PRESSURE  BETWEEN  CHANfaER  and  exit 
C  h=STEP  SIZE  IN  INCHES 

c  csno  of  steps  printed  ie  print  every  l  th  step 
c  MAXsmaXIMUM  no  of  iterations  allowed 
c  n=no  of  steps  taken 

C  IMA=N0  OF  STEPS  PRINTED  -  IE  NA  STEPS 

c  properties  of  tp-hiou  propellant 

C  U=VISCOSlTY  LBF  SEC/FT**2 

C  CiBURNlNG  RATE  COEFFICIENT  FT**(2XP+i)/SEC  LBF**xP 
C  u=PROPELLANT  DENSITY  LbM/FT*w3 
C  10=BURNING  TEMPERATURE  K 
C  GM=RATlO  OF  SPECIFIC  HEATS 
C  aP=bURNING  rate  exponent 
C  h=GAS  CONSTANT  FT  LtiF/LUM  R 
C********t*************************#*****************%*****c 
Cr,0U47 
TOS6200, 

Drill# 

XPS.25 

GMSl.l 

Ur, 00000202 

Rs(l545#)/(26.28) 

C  INPUT  geometry  for  debond 

S(l)rO# 

YP ( i ) r ( YN-Yo ) /LL 
THETA=ATAN(YP(1)) 

Y(l)rYO 

333  DO  2  1 1=2# N 
YP(II)=YP(l» 

S(II)=H*(II-1) 


oooooooooooooooooooo 
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Y ( I I ) sYP ( I I ) *S ( I I ) +Y0 
2  CONTINUE 
Y(N)=YN 

C  print  heading  and  debond  geometry 

998  PRINT  1  LL»Y(1)»Y<N),PC 

999  Nisi 

C  INITIAL  PO  6UESS 
POSPC+1 
PORSsPC 
SRSO. 

60  TO  2b 
24  CONTINUE 
DP0SsP0-P0RS 
DSSS(I)-SR 
IF (DS) 100 » lOlr 100 
101  DS=H 
100  CONTINUE 
SRSS(I) 

PORSsPO 

IF(NI“6) 54»54r53 
54  SDXVS10. 

GO  TO  57 

53  IF  <NI-10.)929»929*930 

929  SD1VS3.0 
GO  TO  57 

930  SDIYS1.0 

57  CONTINUE 

C  CORRECTED  PO  GUESS 

PO=PORS* ( OPOS/SOI V ) * ( S  <  N ) -SR  > /OS 
26  CONTINUE 
BsO. 

AA=(P0*144O**(XP-l.) 

THOsD#C*AA 

T(l)sTHO 

2(l)sl. 

Fs(2.*(GM-l.)*T0*R)/(GM*32.2) 

NlSNl+i 

IF  t  Nl-MAX ) 56  r  58  r  29 

58  CONTINUE 

DO  36  Isl,N 
P(I)sPO  *  2(1) 

IF(I-N)50r22r22 
50  CONTINUE 

IF ( P ( I ) ♦ERROK-PC ) 24  r  22  r  22 
22  CONTINUE 


X ( I )  =SQRT ( 2 ( I ) *Z ( I )  +F#T  ( I ) *T ( I ) ) 
V(I)=<X(I)-Z<I>)/<2(I)*(GM-1.>> 

IFU-N)51, 23*23 
bl  CONTINUE 

IF<V<I)-1.)23»24*2*> 

23  CONTINUE 

RE(I)=(Y(I)*T(I)*P0*24.)/{U*32.2, 

IF  (RE<I>  -  2200.)  20*  20*  21 

FRICTION  FACTORS  FC  ESTIMATED  FROM  REFERENCE  8 

20  FC=64./RE(I) 

60  TO  27 

21  IF  (KE(I)-‘t000.)  400*400*401 

400  FC=.075 
60  TO  27 

401  FC“.07 
27  CONTINUE 

BR  U ) =THO*Z ( I ) **XP*SGRT ( RE  < I >  > /T 1 1 > 

52  CONTINUE 

IF (N-I ) 108*108*105 
105  CONTINUE 

HBH=l./<2.*(bR(I)+l.)m./<2.*5***<BR<J)'4.M 

d2=F*T(I)/X(I) 

b4=(GM*XtI>)/<(X(I)-6M*Z(I)>*Y(I)> 

Bls(X(I)-Z<I))*(FC*(l.+HBR)/2.+YP(I)) 

B3=TH0*Z ( I i **XP-YP (I)*T(I) 

ZP=d4*(Bl+B2*B3) 

2<I+l)=Z(I>+h*ZP 
IF(Z(I+1)+Z(I))24,24*25 
25  CONTINUE  • 

b=(Z(I+i>+2(i))**XP+B 
T (l+l)s(THO/Y(I+l) )*(H*B/2»**XP+Y0> 

108  CONTINUE 

V(I)=SQRT(V(I)) 

36  CONTINUE 

1F(V(N)-0ZF-1. >60*60*59 
60  IF <P(N)+fcRROR-PC >90*56*91 

91  IF ( ABS  <  P ( N ) -PC ) -ERROR  >56*56*92 

92  IF(ABS(VIN>-1.)-DIF)56*56»93 

93  IFIABS<P(N>-PC)/PC-ABS<1.-V(N>> >90*90*88 
59  UPO=PO/100* 

60  TO  89 

90  0P0=(PC-P<N))/10. 

Gg  TO  89 
88  UPO=PO/100» 

69  CONTINUE 


PO=PO+OPO 

EXOPsOPO/lO, 

NV=1 

NP=1 

55  CONTINUE 
VR=V(I) 

PR=P { I ) 

NI=ND1 

IF(Nl-MAX)86»86r29 
86  CONTINUE 
B=0. 

THO=D*C« ( P0*144 • ) ** < XP-1 • ) 

T(l)=THO 

2(1)=1. 

00  37  1=1  »N 

F=(2.*(GM-1.)*T0*R)/<GM*32.2) 

P(I)=PO*Z(I) 

IF(I-N)61*62»62 

61  CONTINUE 

IFCPi I) +ERKOR-PC ) 24 » 62#  62 

62  CONTINUE 

XU)=SQRT<ZU)*Z(I>+F*T<I)*T(n> 

V<I)  =  <X(I)-Z(I>>/<Z(D*(6M-1.)> 

IF(I-N)63, 64*64 

63  CONTINUE 
IF<V(I)-1. >64*64*24 

64  CONTINUE 

RE ( D  =  ( Y ( I ) *  T ( I) *P0*24 • ) / ( U#32*  2 ) 
lF(RE(I)-2200. >65*65*66 

65  FC=64./RE(I) 

GO  TO  67 

66  IF  (RE(I)-4000* )500» 500*501 

500  FC=*075 
GO  TO  67 

501  FCS* 07 

67  CONTINUE 

6R  U  > =TH0*2 ( I ) **XP*SQRT ( RE (I > > /T ( I ) 
HBR=l./(2.*(8R(D+l.))+l./(2.*5.**<BR(D/4.)) 
IF(N-I ) 106# 106* 107 
107  CONTINUE 

62=  F  *  T(I)  /  X(I> 

64  s  (GM*  X(l) )  /  ((X(I)  -  GM  *  Z(I))*Y(I)> 
6l=(X(D-ZtD)*(FC*(l*+HBR)/2.+YP(I)) 

B3=THQ*Z ( I ) **XP-YP ( I ) *T ( I ) 

ZP=B4*(B1+B2*B3) 


£ ( I +1 ) =2 ( I ) +H*2P 
1F(2( I+l)+2( I)  )24#69#69 
69  CONTINUE 

b= ( 2 ( I +1 ) +2 ( 1 ) )  **XP+B 
T(I+l)=(rHO/y(I-H) )*(H*B/2.**XP+Y0) 

106  CONTINUE 

V(I)=S0KT(V(1)) 

37  CONTINUE 
DPsP(N) -PR 
OV=V(N)-VR 

OHC=AbS ( P ( N ) -PC ) /PC-ABS (1 . - V ( N ) ) 
IFtP(N)+ERROH-PC)80,8l»8l 

81  IF(V(N)-LIF-1.)82»56»83 

82  IF(A8S(P(N)-PC)-ERRoR)56»56#87 
87  IF(ABS(V(N)-1.)-DIF)56»56»84 

84  IF(ABS(P(N)-PC)/P(N)-ABS(1.-V<N) ) >80*80*83 
80  CONTINUE 

NP=NP+1 

CNP=ABS(bP/EXOP) 

I F ( NP-2 ) 95 » 9b  »  94 

95  PbIV=3. 

00  TO  99 

94  PDIV=CNP 
99  CONTINUE 

POKzPO- ( IP(N)-PC)/P0XV)*0P0/DP 

EXUP=(P(N)-PC)/PDIV 

GO  TO  85 

83  CONTINUE 
NVSNV+1 

IF ( NV-2 ) 97 » 97  *  96 

96  IF(AtS(V(N)-l* )-«3) 104*103*103 
104  V0IV=2, 

GO  TO  98 

97  IF(ABS(VtN>-l.) -.3)102*103*103 

102  V0IV=3, 

GO  TO  96 

103  VDIV=10.0 

98  CONTINUE 

POHsPO-  (l  V  <  N )  -1 . )  /VI)  I V  >  *DP0/O  V 

85  CONTINUE 
DPO=  POR-PO 
PO=POR 

GO  TO  55 
56  CONTINUE 
601  PRINT  7 


PRINT  8  (S(I),Y(I)»P(I>,V<I>»I=i*NA»L) 

PRINT  PROPERTIES  AT  LOCATION  OP  TRANSDUCERS 
IE  .19  1.72  AND  3.72  INCHES  FROM  TIP 


JJ(l)-20. 

JJ(2)=vJJll>+153. 

J J ( 3 ) s  J  J l 2 ) +200 . 

DO  300  KKsl,3 
KsJJ(KK) 

PRINT  8  S(K),Y(K)*P(K>»V(K) 

300  CONTINUE 

THETA=THETA*57.3 

print  301  xp,r,  u,c,d,gm»theta 


GO  TO  5 
29  CONTINUE 

PRINT  72  w 

1  FORMAT  (lHl,54X,* DEbOND  GEOMETRY* »/» 18X,  Ut|vpu# 

♦  tDEBOND  LENGTH', 13X,*TIP  WIDTH*, 13X, ’EXIT  WIDTH*, 

*  i3x, * chamber  pressure* »/»21x>F5.2,2Qx,F6.4,15X, 

♦  F6.4,20X,F6.1,/») 

format  (28X?*S(I)  inches*, 9X,*YU)  INCHES*, 11X, 

*  »P(I)  PSI*,13X*MACH  NO.*) 

8  FORMAT  (2aX»F10#4,10X,Fl0.4»10X,Fl0.2,10X»F10.4) 

72  FORMATS  *,* ITERATION  UNSTABLE*) 

301  FORMAT  ilH0»10X,'XPs«»F4.2,3X»»Rs»,F4.1»3X,'Us_, 


3 

7 


* 

* 


*vn  r  n*  —  ,  ■  ■  *  _ 

F10.  9,3X,*Cs*,F6.4,3X, *D=» »F5.1,3X» 'GM=» ,F3.1, 

3X, * CRACK  angles* ,F4.2, IX, » DEGREES* ) 


END 


